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SPATIAL STUDY WITH THE VLT 
OF A NEW RESOLVED EDGE-ON CIRCUMSTELLAR DUST DISK 
DISCOVERED AT THE PERIPHERY OF THE p OPHIUCHI DARK CLOUDS 
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ABSTRACT 

We report the discovery in near-infrared (NIR) with SofI at the New Technology Telescope (NTT) 
of a resolved circumstellar dust disk around a 2MASS source at the periphery of the p Ophiuchi dark 
cloud. We present follow-up observations in J, H, and ATg-band obtained with ISAAC at the Very 
Large Telescope (VLT), under 0''4-seeing conditions, which unveil a dark dust lane oriented East-West 
between two characteristic northern and southern reflection nebulae. This new circumstellar dust disk 
has a radius of 2"15 (300 AU at 140pc), and a width of 1'.'2 (170 AU at 140pc). Thanks to its location 
at the periphery of the dense cores, it suffers small foreground visual extinction (Ay = 2.1 ± 2.6mag). 
Although this disk is seen close to edge-on, the two reflection nebulae display very different colors. We 
introduce a new NIR data visualization called “Pixel NIR Color Mapping” (PICMap for short), which 
allows to visualize directly the NIR colors of the nebula pixels. Thanks to this method we identify a 
ridge, (X'S (40 AU at 140 pc) to the north of the dark lane and parallel to it, which displays a NIR color 
excess. This ridge corresponds to an unusual increase of brightness from J to ATs, which is also visible 
in the NTT observation obtained 130 days before the VLT one. We also find that the northern nebula 
shows ~3mag more extinction than the southern nebula. We compute axisymmetric disk models to 
reproduce the VLT scattered light images and the spectral energy distribution (SED) from optical to 
NIR. Our best model, with a disk inclination i = 86 ± 1°, correctly reproduces the extension of the 
southern reflection nebula, but it is not able to reproduce either the observed NIR color excess in the 
northern nebula or the extinction difference between the two reflection nebulae. We discuss the possible 
origin of the peculiar asymmetrical NIR color properties of this object. 

Subject headings: Open clusters and association: individual: p Ophiuchi dark cloud —infrared: stars 
—Stars: individual: 2MASSI1628137-243139 —Stars: pre-main sequence —Stars: 
circumstellar matter —Stars: formation 


1. INTRODUCTION 

Observations of resolved circumstellar disks around 
young stars are very important for our understanding of 
the formation of solar-type stars and planetary systems. 
Usually the large difference of brightness between the cen¬ 
tral star and its circumstellar material does not allow di¬ 
rect imaging of the circumstellar disk. However when the 
disk midplane is seen close to edge-on, models predict 
that the direct light of the star is blocked out while a 
small amount of the light is scattered by the grains be¬ 
low and above the disk midplane, producing a typical pat¬ 
tern of a dark lane between two reflection nebulae (Whit¬ 
ney & Hartmann 1992). The first young stellar object 
displaying this typical pattern, HH30IRS in the Taurus 
dark cloud, was imaged by the HST/WFPC2 (Burrows 
et al. 1996). Edge-on disks have since become accessible 
to ground-based NIR observations thanks to adaptive op¬ 
tics or speckle imaging (e.g., observations of HKTau/c; 
Stapelfeldt et al. 1998; Koresko 1998). Recently new ob¬ 
servations of edge-on disks were reported in one of the 


nearest star-forming regions, the p Ophiuchi dark cloud 
(L1688; d ^ 140pc), both by the VLT/ISAAC under 
O'.'35-seeing conditions (Brandner et al. 2000) and the 
HST/NICMOS (Terebey 2001; Allen et al. 2002). These 
disks are all located in the dense cores of this dark cloud 
and are deeply embedded (Fig. 1). 

During the period April 4-7 2001, five deep NIR point¬ 
ings were obtained with the NTT/SofI in the p Ophi¬ 
uchi star-forming region (Grosso et al., in preparation), 
as follow-up of X-ray observations made with Chandra 
and XMM-Newton. Grosso et al. (2001) already reported 
the serendipitous discovery of new embedded Herbig-Haro 
objects in one of these pointings. We focus here on an¬ 
other pointing at the periphery of the p Ophiuchi dark 
cloud (Fig. 1) where we have discovered a new resolved 
edge-on circumstellar dust disk. We present in §2 the 
NTT/Sofl data and the follow-up observation made with 
the VLT/ISAAC (Unit Telescope 1; Antu) unveiling the 
shape of this object. We measure its photometry in §3, 
and discuss the foreground visual extinction in §4. We 
introduce in §5 the Pixel NIR Color Mapping of the two 
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nebulae. In §6 axisymmetric models are used to estimate 
the physical parameters of this disk, and to predict NIR 
colors. Finally in §7 we discuss the possible origin of the 
peculiar NIR color properties of this object. 
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Fig. 1.— Finding chart of the edge-on circumstellar disks in the p 
Ophiuchi dark cloud. The optical background image from the Dig¬ 
itized Sky Survey is merged with the ISOCAM mid-infrared map 
(6.7 and 14.3 /im; Abergel et al. 1996). The overlaid white contours 
show the column density (Av=36, 54, 72, 90mag; Wilking 

Sz. Lada 1983), associated with the L1688 dark cloud. Black dia¬ 
monds mark the previously known young stellar objects displaying 
in scattered light two reflection nebulae separated by a dark lane 
(Brandner et al. 2000; Terebey 2001; Allen et al. 2002). They are 
mostly protostars surrounded by remnant infalling envelopes. The 
dashed box displays the 5^x5^ field of view of the NTT/SofI point¬ 
ing at the periphery of the dark cloud. The small square box defines 
the area of Fig. 2. 
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2. OBSERVATIONS 
2.1. The NTT/SofI observation 

This observation was made during the night of April 
7 at airmass 1.0, with the large SofI field (5' x 5'; pixel 
size = 07292) to survey the largest area. We took 60 s- 
exposure frames using 10 s detector sub-integration time 
to avoid detector non-linearity and saturation. Between 
each frames the telescope was moved inside a 20"-width 
box centered on the target position, according to a ran¬ 
dom pattern of offsets automatically determined by the 
imaging mode called auto jitter. The net exposure times 
are 24min in J, H, and Kg, leading to a ATs^lSmag 
completeness magnitude for point sources. The data re¬ 
duction was performed using the ESO’s eclipse pack¬ 
age (version 3.8.1) which provides the jitter command. 
The best seeing conditions of this NTT/SofI pointing were 
achieved during the Ks exposure with Full Width at Half 
Maximum (FWHM) ~(y.'7. The astrometric registration 
of the reduced images was made using the position of a 
773=14.1 mag star (2MASSI1628133-243127; Cutri et al. 
2000) located only 12" to the north of our object, extracted 
using SExtractor (Bertin & Arnouts 1996). J and H im¬ 
ages were then resampled to match the 77s image using the 
warping command of eclipse with the hyperbolic tangent 
interpolation kernel. 

Figure 2.a shows the resulting color composite image 
centered on 2MASSI1628137-243139 which went unno¬ 
ticed until now. Compared to the unresolved field stars, 
this source appears extended. One can barely see two 
characteristic reflection nebulae, and a marginally resolved 
^ http://www.eso.org/outreach/press-rel/pr-2002/pr-09-02.html 


dark dust lane oriented East-West in front of the star. To 
confirm this detection and to measure the physical pa¬ 
rameters of this disk, we need a deeper and higher angular 
resolution image. This can easily be achieved with the 
VLT-UTl/ISAAC, by combining better seeing conditions 
with better sensitivity and image sampling. We obtained 
for this goal one hour of the ESO Director’s Discretionary 
Time. 


2.2. The VLT-UTl/ISAAC observation 

The observation was carried out in service mode, trig¬ 
gered by very good seeing conditions (FWHM~074) on 
August 15, 2001. The auto jitter imaging mode was also 
used with 80s, 60s, 60s-exposure frames and 40s, 15s, 
10s detector sub-integration time for J, 77, and 77s, re¬ 
spectively. Net exposure times were 10 min in Kg and 77, 
and 10.7min in J, with airmass ranges from 1.05 to 1.11. 
The data reduction was performed using the version 4.1.0 
of eclipse. We used the is_gliost command to remove 
the electrical ghosts induced by a bright star in the field of 
view. The photometric calibration is based on the obser¬ 
vation of the IR standard star 9181/S243-E (Persson et al. 
1998) provided by the service mode, and taking into ac¬ 
count average IR atmospheric extinction on Paranal (0.11, 
0.06, and 0.07 mag airmass”^ for J, 77, and 77s, respec¬ 
tively). We used the same method as described above for 
our NTT data to produce the resulting color composite im¬ 
age presented in Fig. 2.b (see also the ESO Press Release, 
Photo 12b/02^, which presents a larger color composite 
image corresponding to the field of view shown in Fig. 5). 

Thanks to very good seeing conditions and smaller de¬ 
tector pixel size, the two reflection nebulae are unveiled, 
and the dark dust lane is now clearly resolved in all three 
filters. The shape of this object seen nearly edge-on and 
appearing on a starry background has led us to nickname 
it the Flying Saucer (hereafter the FS). Two sources are 
visible within 2" (280 AU at 140 pc) of this disk both in 
the NTT and the VLT images. The first source located 
at the North-East is the brightest with Kg ~ 17.5 mag. 
Stars with similar colors are visible in Fig. 2.b; this source 
is likely a background star. The second source, thanks to 
its location close to the South-North disk axis, could be 
a jet H 2 knot candidate, although no counter knot is vis¬ 
ible. Several objects with similar colors and brightnesses 
are also spread over the field of view, which allows us to 
interpret this second source also as a background source. 

3. PHOTOMETRY OF THE FS 

Figure 3 presents the surface brightness maps of the 
VLT data in the J, 77, and 77s-band, respectively, with 
contours obtained without any smoothing with the Inter¬ 
active Data Language (IDL). We define the dark lane axis 
as the bisector of the line segment defined by the nebula 
peaks in the Kg image. The orientation of the dark lane is 
strictly East-West, and separates the northern and south¬ 
ern reflection nebulae. The southern nebula is slightly 
asymmetric, with a larger extension towards the East. In 
J the southern nebula is more extended than its northern 
counterpart, and the southern peak is 1.7 times brighter 
than the northern one. This implies that we are not seeing 
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Fig. 2.— Near-infrared view of 2MASSI1628137-243139. These color composite images combine J, and using respectively blue, 
green, and red color coding, with logarithmic color stretch. The white balance is obtained from the brightest background star. The image 
size is 30^^x30^^ The image orientation is North on the top and East on the left hand side, a) observations with the NTT/Sofl (pixel size = 
0'.'292, seeing FWHM~0'.'7). b) observations with the VLT-UTl/lSAAC (pixel size = 0'.'148, seeing FWHM~0'.'4). 


the FS exactly edge-on. Indeed the scattered light models 
of accretion disks versus inclination show that when one 
looks above the disk midplane, the brightness of the up¬ 
per reflection nebula increases compared to the brightness 
of the lower reflection nebula (see Whitney & Hartmann 
1992). We are therefore seeing the FS slightly below its 
disk midplane, corresponding to the southern nebula. The 
difference between the northern and the southern reflec¬ 
tion nebula decreases with increasing wavelength, leading 
to a more symmetric reflection nebulae in K^. In this fil¬ 
ter the northern nebula is still slightly less extended, but 
its peak brightness is 1.3 times brighter than the southern 
nebula. This produces the reddish aspect of the northern 
nebula visible in Fig. 2 both in the NTT and the VLT 
images. This resolved circumstellar disk has a radius of 
2.15" (300 AU at 140 pc), as measured from the largest ex¬ 
tension of the southern nebula (contour 20.4 mag arcsec“^ 
in J), and a width of 1.2" (170 AU at 140pc). The shape 
of the FS, due to its low disk inclination, is reminiscent 
of HKTau/c (Stapelfeldt et al. 1998; Koresko 1998), but 
three times larger and without strong asymmetry in the 
scattered light intensity. 

We compute the photometry of the integrated light of 
the FS using an aperture identical for the three filters, 
and background estimate from an area free of sources 
close to the disk. To avoid contamination by the close 
background star, the aperture photometry is computed 
from the area where the surface brightness is greater than 
19.3 mag arcsec“^ in the three filters (see the boundary 
of the aperture photometry in Fig. 6). The match be¬ 
tween the ISAAC filters used here and those used to es¬ 
tablish the faint IR standard star system of Persson et 
al. (1998; a.k.a. the LCO photometric system) being quite 


good, we use the color transformations given by Carpen¬ 
ter (2001) for the LCO photometric system, to convert the 
disk magnitudes to the 2MASS photometric system. The 
resulting J, H, and Ks integrated magnitudes are listed 
in Table 1. These values are consistent with the values 
given for 2MASSI1628137-243139 in Cutri et al. (2000), 
taking into account the photometry uncertainties. This 
implies no dramatic change in the brightness of this disk 
on ~3 years timescale. However this does not exclude very 
large changes over a few days like those seen in HH 30 IRS 
(Wood et al. 2002). 

We list also in Table 1 existing complementary inte¬ 
grated magnitudes from other broad-band optical/NIR Al¬ 
ters. Using the absolute calibration given by Cox (2000), 
we obtain the corresponding spectral energy distribution 
(SED) from optical to NIR presented in Fig. 3. This disk 
was not detected in the mid-infrared (MIR) by the ISO¬ 
CAM survey although located in low-noise region of this 
survey (Bontemps et al. 2001; 6 and 10.5 mJy upper lim¬ 
its for 6.7 and 14.3/im, respectively, corresponding to the 
central wavelength of the ISOCAM filters LW2 and LW3). 
Indeed for edge-on inclinations, all the material of the disk 
midplane is on the line of sight leading to an optically thick 
medium in the MIR; continuum emission from the warm 
inner disk needs deeper exposure times to be detected (see 
e.g., ISOCAM MIR spectrum of HH30IRS; Stapelfeldt & 
Moneti 1999). The SED shape, rising in the optical and 
reaching a peak in H, is similar to that of HH 30 IRS (opti¬ 
cal variability from Wood et al. 2000; NIR from Cotera et 
al. 2001; mid-infrared to millimetric from Brandner et al. 
2000 and references therein), and consistent with edge-on 
disk models (e.g. Boss & Yorke 1996; see also below §6.4). 
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Fig. 3. — Surface brightness maps of the Flying Saucer from the VLT-UTl/ISAAC observations. We give the calibrated contour intensity 
range from the peak value to the first contour, with the contour step. The crosses mark the peak positions of the reflection nebulae in each 
filter. The symmetry axes (dashed lines) are defined from the peak positions in the JFs-band image. Their lengths give the width and the 
height of the disk. 


Table 1 

Optical/NIR integrated photometry of the Flying Sauger. 


a 

b 

a 

C 

d 





B 

-^53F 

R 

I 

J 

H 


20.4±0.3 

20.5±0.3 

18.6±0.3 

18.4±0.3 

17.16±0.14 

15.93±0.06 

14.78±0.05 

13.98±0.03 


®SuperCOSMOS Sky Survey (Hambly et al. 2001). 

^Derived from Bj (Kodak Illa-J plate + GG395 filter) using the color correction of Blair & Gilmore (1982) assuming the color index V — R 
corresponding to = 3500 K in the MK temperature scale (Cox 2000). This leads to V = 19.9 (marked with a white dot in Fig. 3). 
‘^Derived from i? 53 F (Kodak IIIa-F plate + RG630 filter) using the color corrections of Bessel (1986) with I given by DENIS (Epchtein et al. 
1999). 

‘^DENIS photometry (M. H. Vuong, private communication). 
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Fig. 4.— Spectral energy distributions (SED) of the Flying Saucer 
compared to others circumstellar disks. The continuous line and the 
dots show the SED of the Flying Saucer with two-sigma error bars. 
(The white dot is only an estimate of the K-band magnitude, see 
note b of Table 1.) The SED of CRBR 2422.8-3423 located in one 
of the dense core of the p Ophiuchi cloud (Brandner et al. 2000), 
and HH 30 IRS located in the Taurus dark cloud, are plotted for 
comparison. 

Since the FS falls outside the region mapped in the mil- 
limetric continuum by Motte, Andre, & Neri (1998) we 


cannot have a direct measurement of the amount of cir¬ 
cumstellar material. An estimate of the disk mass will be 
given from scattered light model in §6. 

4. FOREGROUND VISUAL EXTINCTION 

To estimate the foreground visual extinction of our ob¬ 
ject, Ay, we have to deredden its integrated NIR col¬ 
ors, derived from the integrated magnitudes established 
in the previous section. Meyer, Calvet, & Hillenbrand 
(1997) observed in Taurus-Auriga that dereddened Clas¬ 
sical T Tauri stars (CTTS) follow a narrow segment locus 
in a J — H,H — K color-color diagram. We adopt this lo¬ 
cus as statistically representative of the CTTS colors. We 
stress that an estimate of Ay based on these typical CTTS 
colors is however only indicative because the NIR light of 
the FS is almost purely scattered. Its color might then 
be quite different from that of a typical CTTS in which 
one primarily sees direct light from the stellar photosphere 
and the inner disk. The color transformation of the CTTS 
color locus from the CIT to the 2MASS photometric sys¬ 
tem using Carpenter (2001) leads to: 

(J - H)o = 0.61 ± 0.12 (iJ - Ks)o + 0.50 ± 0.07. (1) 

Meyer, Calvet, & Hillenbrand (1997) used the extinction 
law from Cohen et al. (1981), without taking into account 
its uncertainties, to deredden their CTTS sample. To be 
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Fig. 5. — Extinction map of the area around the Flying Saucer. The background image is the Xs-band VLT-UTl/ISAAC image displayed 
using a logarithmic color stretch. The dashed box corresponds to the area of Fig. 2. Squares and ellipses mark stars and galaxies, respectively, 
with signal-to-noise greater than 5, selected using SExtractor (Bertin & Arnouts 1996). Contours show the extinction map computed from 
the colors of the selected (background) stars. Contour labels give the visual extinction in magnitude. 


consistent with their results we take the same extinction 
law transformed to the 2MASS photometric system: 

(J-iJ) - (J-iJ)o = 0.118 ±0.001 Av, (2) 
{H - Ks) -{H - Ks)o = 0.067 ± 0.001 Ay, (3) 

which leads to E{J - H)/E{H - Ks) = 1.77 ± 0.04. 
Straightforward algebra with formulae 1-3 using the ob¬ 
served NIR colors {J — El = 1.15 ± 0.08 mag, H — Ks = 
0.80 ± 0.06mag), and combining the quoted uncertainties, 
gives 2lv=2.1±2.6mag (see Fig. 6). 

In order to obtain an upper limit on the foreground 
visual extinction of the FS, we computed the extinction 
map of the area around it by applying the “Near Infrared 
Color-Excess” (NICE) method (Lada et al. 1994; Alves 
et al. 1998). We extracted with SExtractor (Bertin & 
Arnouts 1996) sources with signal-to-noise greater than 5 
in the three filters. To compute this map we selected only 
sources with stellarity indexes greater than 0.95 (other 
sources were considered as galaxies). Because there were 
no observations of the unreddened background stellar field 
(necessary for the determination of the zero point of the 
extinction scale), we must adopt a < H — Ks > color for 
this background. Alves et al. (1998) found < H — Ks > 


= 0.20 ± 0.13 mag towards a complex region in Cygnus, 
essentially at the Galactic plane, while Alves, Lada, & 
Lada (2001) found < H — Ks >= 0.12 ± 0.08mag to¬ 
wards a clean line-of-sight towards the Galactic bulge, at 
5 ~ 7°. The stellar background towards the Ophiuchus 
complex is probably even better behaved, because it lies 
at 5 ~ 15°. On the other hand, the < H—Ks > color of the 
North Galactic pole from 2MASS data is also ~0.12mag 
(M. Lombardi, private communication), so we decided to 
adopt this as the mean Ophiuchus background color. 

We show our extinction map in Eig. 5. The spatial 
resolution of the map is 20'', and the visual extinction 
ranges from 5 to 7 mag. The extinction increases from 
East to West, i.e. towards the molecular core (see Eig. 1), 
hence tracing the presence of material related with the 
main cloud. This extinction map leads to the follow¬ 
ing upper limit on the foreground extinction of the FS: 
Ay ± 5.1 ± 1.2 mag. Our previous estimate based on the 
CTTS colors, Av=2.1±2.6mag, is consistent with this up¬ 
per limit, and might imply that the FS is located within 
this peripherical material, which is not unreasonable. We 
adopt this extinction value as our best guess. 
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Thanks to its location at the periphery of the dense 
cores, this disk suffers less foreground extinction than the 
other disks already found in this star-forming region. This 
is also visible in Fig. 3 where its SED shows only a slow 
rise at short wavelength in contrast to the steep rise of 
CRBR 2422.8-3423 due to the high extinction of dense 
cores (Brandner et al. 2000). Cotera et al. (2001) found 
a comparable foreground visual extinction for HH 30 IRS 
{Ay < 4.1 ± 0.9 mag) assuming a typical stellar type MO 
for the central star. Since the FS suffers no background 
emission from the dense cores, it is an ideal laboratory for 
the study of both dust and gas in circumstellar disks. 

5. PIXEL NIR COLOR MAPPING 

OF THE REFLECTION NEBULAE 

We noted in the previous section the reddish appear¬ 
ance of the northern reflexion nebula of the FS, appearing 
in Fig. 2. To our knowledge such a NIR color difference 
between the reflexion nebulae of an edge-on disk has never 
been observed. NIR adaptive optics images of HK Tau/c 
display only a “slight change in the nebula’s appearance 
between J and K” (Stapelfeldt et al. 1998); speckle imag¬ 
ing shows that the flux ratio between the two nebulae in¬ 
creases from 5 at J to 8 at if (Koresko 1998). To do a 
spatially resolved quantitative study of the NIR color dif¬ 
ferences between the two nebulae of the FS, we introduce 
a new NIR data visualization method called “Pixel NIR 
color mapping”, PICMap for short. Although in this paper 
we apply it to the NIR photometry of the FS, this method 
is general and could conceivably be applied to other pho¬ 
tometric bands and adapted to other spatially resolved 
objects. 

The central idea is the following. We want to visualize 
directly in a single VLT image the J — H and H — Kg 
colors of each pixel of the FS. We start by selecting the 
pixels of the FS located inside the boundary of the aper¬ 
ture photometry defined in §3 (see the upper left-panel 
of Fig. 6), and by computing their J — H and H — Kg 
colors from the J,H and Kg brightness maps of Fig. 3. 
Then we plot this set oi J — H and H — Kg colors in a 
J — H, H — Kg color-color diagram, much in the way the 
NIR colors of individual stars in an association would be 
plotted. By analogy, we also plot in this diagram the stan¬ 
dard ISM reddening vector. The upper right-hand panel 
of Fig. 6 illustrates the main result: most of the pixels 
in the J — H,H — Kg diagram are found to be spread 
roughly along the reddening vector. Therefore, we use its 
direction to build an arbitrary color palette, shown to the 
right of the diagram. The last step of the method is then 
to create a mapping, both on the J — H, H — Ks diagram 
and on the VLT image, by “coding” the color of each FS 
pixel according to its projection on this color palette. In 
other words, PICMap defines a correspondence between 
the J — H, H — Ks plane and the VLT image. As we shall 
see, this visualization method is powerful to detect any 
unusual trend in NIR colors within the FS image. 

Physically, we could interpret the pixel color palette 
used in the PICMap as a visualization of the NIR col¬ 
ors of the nebula pixels in terms of pure extinction by 
interstellar medium type dust grains. We caution however 
that the grains in the disk may be grayer than interstellar 
medium type dust grains (see §6.2), and the true extinc¬ 


tion along a line-of-sight through the nebulae is then likely 
considerably higher than the standard extinction derived 
from this simple ISM extinction law. In reality, the NIR 
colors of the FS result from a complex combination of ex¬ 
tinction and scattering by disk grains, which can only be 
reproduced by modelling the radiative transfer of the stel¬ 
lar light through the dust, as discussed in the next section. 
However the PICMap can be used to quantify the relative 
difference of NIR colors between the two reflection nebulae 
in term of standard ISM extinction. 

In the upper panels of Fig. 6 when one moves in the 
VLT image from the edge of the southern nebula to the 
dark lane along the vertical axis of the disk, the corre¬ 
sponding NIR colors in the J — H, H — Kg diagram (white 
dots) follow roughly the direction of the ISM reddening 
vector, with an increase of ~4mag of standard extinction. 
The same amount of standard extinction is also visible for 
the NIR colors of the pixels located between the edge of 
the northern nebula and the dark lane (black dots). How¬ 
ever the northern nebula suffers an extra standard extinc¬ 
tion of ~3mag compared to the southern nebula. In the 
J — H, H — Kg diagram, it is also clear that the NIR colors 
are at a peak for many pixels of the northern nebula, and 
display what would be called a “NIR color excess” in the 
case of young stars. The PICMap method allows to locate 
them easily spatially on the VLT image. They are not 
distributed at random in the northern nebula, but form 
a ridge, 0'.'3 (40 AU at 140 pc) to the north of the dark 
lane and parallel to it. There is no similar feature in the 
southern nebula. 

In spite of the fact that the FS disk is seen nearly edge- 
on, it is clearly asymmetric with respect to the dark lane, 
the northern nebula being significantly redder than the 
southern nebula. In the next section, we will use axisym- 
metric disk models with radiative transfer, to look for a 
model parameter set that reproduces the observed NIR 
color asymmetry. We will apply our PICMap method to 
the scattered light model, to visualize the results and com¬ 
pare it with our observations. 

6. RADIATIVE TRANSFER MODELS 

We adopt a axisymmetric density structure for the cir¬ 
cumstellar material (i.e. symmetric with the disk mid¬ 
plane), with dust grains that reproduce the images and 
spectral energy distribution (SED) of HH30 IRS, and per¬ 
form the radiative transfer simulation using Monte Carlo 
techniques (e.g., Whitney & Hartmann 1992; Code & 
Whitney 1995) to compare scattered light models with 
the high resolution VLT images. For calculating synthetic 
SED, we use the Monte Carlo radiative equilibrium tech¬ 
nique of Bjorkman & Wood (2001) as adapted for T Tauri 
disks (Wood et al. 2002a, b). 

6.1. Density structure of the circumstellar material 
6.1.1. Disk model 

As with previous models of disk SEDs and scattered 
light images, we adopt a axisymmetric flared disk struc¬ 
ture in which the total density (gas plus dust) is: 

«,.. = «(!-) exp (-^^), (4) 

where po is the disk density extrapolated to the stellar 
surface (we set i?* = 1.03 R©; see SED model §6.4), 
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Fig. 6.— Pixel NIR Color Mapping. Upper panels: a) (H — Ks,J — H) map of the VLT images. The background image is an enlargement 
of the color composite image presented in Fig. 2.b, resized to the field of view of Fig. 3. The white crosses mark the peak positions of the 
reflection nebulae in the it's-band image, and the dashed lines are the symmetry axes drawn in Fig. 3. The white contour defines the boundary 
of the aperture photometry defined by surface brightness at least greater than 19.3 mag arcsec”^ in the three NIR filters. The color of each 
pixel of this area is coded according to its position in the J — H, H — Ks diagram on the right-hand panel. The purple arrow shows the path 
when one moves on the vertical axis of the disk; the corresponding path is also shown in the right-hand panel, b) Corresponding J — H, H — Kg 
diagram of the nebula pixels (2MASS photometric system). The black and dotted continuous lines show the intrinsic colors of A0-M6 dwarfs 
and giants (Bessel Sz Brett 1988), respectively. The orange line is the CTTS locus from Meyer et al. (1997). Reddening vectors (Cohen et al. 
1981), adapted for the 2MASS photometric system, are drawn for giants, M6 stars, and CTTS. The orange diamond shows the NIR colors of 
the integrated light of the disk. The arrow marks the estimate of the disk foreground visual extinction obtained from the CTTS locus (see 
§4), and gives the extinction scale. The white and black squares mark the NIR colors of the integrated light of the southern and northern 
nebula, respectively. White and black dots mark the NIR colors of pixels of the southern and northern nebula, respectively. The one-sigma 
photometric errors are given by the sizes of the error boxes centered on the dots. The color palette defined on the right parameterizes the 
NIR colors of the nebula pixels in terms of simple ISM extinction. The color of the error boxes are coded according to their projections on 
the color palette. The purple arrow corresponds to the NIR color variations when one moves on the vertical axis of the disk. See detailed 
discussion of this Fig. in §5. Lower panels: a) [H — Kq,J — H) map of the disk plus diffuse infalling envelope model, b) Corresponding 
J — — Kq diagram of the model image pixels. Symbols are the same than used in the upper panels. The color palette is identical to the 

one use in the upper panel, allowing easier comparisons between the model predictions and the observations. 
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w is the radial coordinate in the disk midplane (with 
w > i?*), and the scale height increases with radius, 
h{w) = ho . We assume that the gas and dust 

are well mixed throughout the disk. We adopt for the flar¬ 
ing parameter (3 the value derived in the hydrostatic disk 
structure models of D’Alessio et al. (1998), (3 = 5/4; and 
set a = 3{P — 0.5) = 9/4, which is the value appropriate 
for viscous accretion theory (Shakura & Sunyaev 1973). 
This leads to a surface density S oc (D’Alessio et al. 
1998), flatter than the one usually found in the literature 
(E (X e.g., Beckwith et al. 1990). The disk outer 

radius, vuout, is fixed to 300 AU to match the extent found 
in the VLT J image. This value is a lower limit of the 
disk outer radius. In our model the disk is passive, with¬ 
out any source of internal heating. For low mass disks 
(Afdisk << 0.1 Mq), the heating due to accretion luminos¬ 
ity is negligible and stellar irradiation dominates the disk 
heating (D’Alessio et al. 1998). The light coming out of 
this system is thus only the light of the star reprocessed 
by the disk dust grains. 


6.1.2. Diffuse infalling envelope model 


It is necessary to add another circumstellar component 
because a model featuring only by a disk produces re¬ 
flection nebulae with a flat hourglass shape, without the 
rounded tips visible in Fig. 3. Inspired by the modelling of 
HH30IRS (Wood et al. 1998), we have included a diffuse 
infalling envelope in our model. The envelope structure 
follows the rotational collapse geometry of Terebey, Shu, 
& Cassen (1984) in which the total density (gas plus dust) 
is: 


Penw 


Menv 

47r 




1 . 

2 



VMo r 


-1 


3 

2 



( 5 ) 


where Menv is the mass infall rate, Vc is the centrifugal 
radius, pL = cosd (with 9 the polar angle measured from 
the disk axis), and pio = P'ir oo) is determined by 


pt+p.o(—-Ij-fi—=0. ( 6 ) 

\ 

For our models we adopt Vc = tUout = 300 AU, M* = 
0.5 Mq, and an outer envelope radius of 1000 AU. We also 
include a curved bipolar evacuated cavity in the envelope. 
The cavity walls are curved with their height scaling as 
-Zcav ~ and a half opening angle of 20° at the outer 
radius of the envelope. 


6.2. Properties of the disk dust grains 

In young stellar objects, the wavelength dependence of 
the dust properties are inferred from the slope of the sub- 
millimetric continuum emission (e.g., Beckwith et al. 1990; 
Beckwith & Sargent 1991) and also from the wavelength 
dependence of the width of dust lanes of edge-on disks 
(e.g., Cotera et al. 2001). These studies show that opac¬ 
ity laws shallower than the one of the ISM are needed. 
Indeed for HH 30 IRS, ISM type grains produce very com¬ 
pact images in it's (Wood et al. 1998), without dark lanes. 
We use the dust grains described in Wood et al. (2002), 


which reproduce the multi-wavelength images and contin¬ 
uum fluxes for the HH 30 IRS disk. This dust model uses 
a “power law with exponential cutoff” for the dust size 
distribution. The exponential scale length for this dust 
distribution, Oc = 50 /rm, indicates dust grains have grown 
to larger than 50 /xm within the HH 30 IRS disk. 

The parameters for the gas plus dust mixture for the 
NIR scattered light models are the opacity, k, albedo, w, 
and Henyey-Greenstein scattering asymmetry parameter, 
g (Henyey & Greenstein 1941), and are shown in Table 2. 
The opacity for this model varies from kj = 24cm^g“^ 
to kk = 16cm^g“^, while for the ISM the opacity is 
larger and has a steeper wavelength dependence with 
Kj = 65cm^ g“^ and kk = 22 cm^ g“^ (e.g-, Kim, Martin, 
& Hendry 1994). 


Table 2 

Parameters adopted 


FOR THE 

GAS PLUS 

DUST 

MIXTURE 

Filter 

band 

a 

K 

[cm^ g“^] 

b 

UJ 

c 

9 

J 

24 

0.53 

0.59 

H 

19 

0.53 

0.58 

K 

16 

0.54 

0.57 


6.3. Scattered light model of the FS 

Having set up the circumstellar material density struc¬ 
ture and dust grain properties as used for modelling other 
circumstellar disks, we compute a synthetic scattered light 
image at A=1.25/xm, convolved with a O'.'4-FWHM Gaus¬ 
sian point spread function. We do not include additional 
noise. We vary the disk inclination, i, the disk scale height, 
ho, disk mass, Mdisk, and envelope accretion rate Menv 
until we reproduce the VLT J image. Then we make scat¬ 
tered light models for the central wavelengths of the H 
and Kg filter. The intrinsic luminosity and spectrum of 
the central star is not set in this calculation. Instead the 
model surface brightness is normalized in each filter us¬ 
ing the peak surface brightness of the southern nebula. 
Among all the free parameters of our model we note that 
only the disk inclination could produce an asymmetry be¬ 
tween the two reflection nebulae in the resulting scattered 
light images. 

Figure 7 shows our simulation for a 300 AU outer ra¬ 
dius disk viewed at z = 86° with ho = 0.015 R*, giving 
h(100AU)=15.3AU, and Mdisk = 21O-3M0. For com¬ 
parison Wood et al. (2002) find using the same modelling a 
similar disk mass for HH 30 IRS with Mdisk = 1-5 10“^ M©. 
From the various models computed with the gas plus 
dust mixture parameters as given in Table 2 we estimate 
the following uncertainties: a-, = 1°, crh(iooAU) = lAU, 
= O.51O“^M0. The disk mass we derive is ob¬ 
viously sensitive to the dust grain properties. For in¬ 
stance for fixed values of w and g, what is really deter¬ 
mined from the modelling is the product of disk mass and 
opacity. If we were to use a different dust opacity the 
disk mass would scale according to that. Moreover the 
disk mass is a lower limit since the opacity is not sensi¬ 
tive to very large particles such as rocks, planetesimals, 
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Fig. 7.— Surface brightness maps of our disk model (see text for the parameter values). 


or planets. The envelope is very tenuous, with a total 
mass Menv = 4 X 10“^ M©. The envelope accretion rate 
is very weak with Menv = 10“^ M© yr“^. For compari¬ 
son typical envelope accretion rates of evolved protostars, 
Class I sources (Lada 1991), are larger with Menv = (2- 
10) X 10“®M©yr“^ (Strom 1994; Kenyon et al. 1993; 
Whitney, Kenyon, & Gomez 1997), whereas CTTS have 
ATdisk = 10“®-10“® M© yr“^ (Strom 1994, and references 
therein). The circumstellar material of the FS is thus char¬ 
acteristic of a CTTS. Our symmetric disk model repro¬ 
duces correctly the North-South extent of the FS in the 
three bands thanks to this diffuse infalling envelope, and 
the East-West shape of the southern nebula. However the 
width of the northern nebula appears ~ 25% larger than 
in the VLT data. A smaller inclination (~ 80°) improves 
the match with the observation by reducing the width of 
the northern nebula, but decreases also dramatically its 
peak brightness leading to a more important discrepancy. 

There are also discrepancies between our axisymmetric 
model and the VLT observation as we go to longer wave¬ 
lengths. In particular the width of the dark lane is repro¬ 
duced only in the J-band. At longer wavelengths the width 
of the dark lane in the axisymmetric models decreases as 
the dust opacity decreases, and the brightness of the north¬ 
ern nebula is always less than that of the southern nebula 
(e.g., Wood et al. 1998). However, the VLT data show no 
variation of the width of the dust lane with the wavelength, 
and moreover as already noted in §3 the northern nebula 
brightness relative to the southern nebula brightness in¬ 
creases. A steeper grain size distribution and/or smaller 
maximum grain size may explain the absence of variation 
of the width of the dark lane at NIR (D’Alessio, Calvet, 
& Hartmann 2001), but cannot reproduce the brightening 
of the northern nebula. The determination of an adapted 
dust opacity law for this object is beyond the scope of this 
paper, and should be inferred using complementary obser¬ 
vation made at longer wavelengths, more constraining for 
the dust size distribution (e.g., Wood et al. 2002). 

The lower panels of Fig. 6 shows the resulting PICMap 
for our model. We can see easily with this tool that 
the NIR color distribution of the southern nebula along 
the reddening vector is quite well reproduced, however 
the NIR color distribution of the two nebulae is nearly 


symmetric with the dark lane. We conclude that the ax¬ 
isymmetric disk model does not reproduce the asymmetry 
in the color spatial distribution observed in the scattered 
light of the FS. 

6.4. BED model of the FS 

Since the spectral type of the central star is unknown, 
we have assumed Teff = 3500 K, and used the correspond¬ 
ing Kurucz atmosphere model (Kurucz 1994) for the input 
stellar spectrum. As the envelope has a small mass com¬ 
pared to the disk mass, distributed throughout a much 
larger volume than the denser disk, we neglect it and 
we compute the synthetic emergent SED from our sym¬ 
metric disk model only. We apply the foreground vi¬ 
sual extinction estimated in §4 using the analytical fits 
from optical to IR of the extinction cross-section given by 
Ryter (1996), combined with the canonical relation be¬ 
tween visual extinction and column density obtained by 
Predehl & Schmitt (1995). Adopting a source distance 
d = 140 pc we require that L* = O.ldL© (which corre¬ 
sponds to i?* = 1.03 R©) to reproduce both the optical 
and the NIR fluxes (Fig. 8). Below 2 fj. the SED is domi¬ 
nated by the (scattered) light of star, while at longer wave¬ 
lengths the emission from the disk dust grains dominates. 
Varying the disk inclination from pole-on to edge-on re¬ 
duces the flux because the disk midplane is optically thick 
in the MIR. Our SED simulation predicts 10-20 /rm fluxes 
of 0.7-1.8mJy, and 1.3 mm fluxes of 12.5 mJy. 
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Fig. 8.— Spectral energy distribution for the Flying Saucer model. 
The black dots mark the observed optical and NIR fluxes with two- 
sigma error bars. (The white dot is only an estimate of the V-band 
magnitude, see note b of Table 1.) The continuous line shows our 
best model: an edge-on disk (short-dashed line) plus 2.1 mag fore¬ 
ground visual extinction. The dotted line is the Kurucz atmosphere 
model for Teff = 3500 K taken as the spectrum of the central star. 
The long-dashed line is the emergent SED when the disk is viewed 
pole-on. 

7. DISCUSSION 

The spatial study of circumstellar disk is usually based 
only on surface brightness maps as in Fig. 3. Thanks to 
the high spatial resolution achieved with our VLT NIR ob¬ 
servations (O'.'4-seeing images), we have been able to make 
a spatial study of the NIR colors of the individual pixels 
of the two resolved reflection nebulae of this nearly edge- 
on disk. We introduced in this paper a new NIR data 
visualization called “Pixel NIR color mapping” (PICMap 
for short), to visualize directly the NIR colors of the neb¬ 
ula pixels. This is a power tool to identify and localize 
any unusual NIR color trend in the reflection nebulae. We 
have found with it a NIR color excess and an extra stan¬ 
dard extinction in the northern reflection nebula. Fitting 
the images and the SED with axisymmetric disk model 
leads to model parameters consistent with what would be 
expected for a typical CTTS, but cannot reproduce the 
observed NIR colors. It is necessary to introduce a source 
of asymmetry in our axisymmetric modelling. 

We use in our model a uniformly bright star. There 
could be cool or dark spots on the star that will change 
the illumination pattern and consequently the brightening 
of the northern nebula. Stellar hotspot models produce 
a typical lighthouse effect, which has been proposed to 
explain the observed brightening and dimming of the re¬ 
flection nebulae in HH30IRS (Wood & Whitney 1998). 
However such a phenomenon should produce changes at 
least on a stellar rotation period (~few days), it is thus 
unlikely to be observed at the same phase on a 130 day 
time interval. We could have an embedded shock associ¬ 
ated with a micro-jet, which could emit H 2 lines producing 
the brightening in H and Kg (e.g., Grosso et al. 2001). 
However there is no evidence in this object of H 2 emission 
knots outside the nebulae. This micro-jet hypothesis in 
the northern nebula can be easily tested with NIR spec¬ 
troscopy and narrow band filter imaging. 

We could account for the lower brightness of the north¬ 
ern nebula at short wavelengths and its redder infrared 
color as due to ~ 3 mag of localized foreground extinction. 


From the color-color diagram of the model image (see lower 
right-hand panel of Fig. 6), we note that adding an extra 
extinction of ~ 3 mag only in front of the northern nebula 
reproduces the observed range oi J — H colors. However 
this localized foreground extinction is not visible in our 
extinction map (Fig. 5). Adding ~ 3 mag foreground ex¬ 
tinction only above the dark lane can only be an ad hoc 
solution. A more natural way to explain these differences 
would be to introduce an extra self-absorption, possibly 
due to a large-scale disk warp. For instance Launhardt 
& Sargent (2001) have detected the millimeter range the 
presence of a symmetric 20° warp beyond 120 AU in an 
~ 200AU-radius edge-on disk. If a such large-scale disk 
warp is seen in front of the central star, it could strongly af¬ 
fect the scattered light. Scattered light models of warped 
disks are now needed to obtain quantitative predictions. 
This idea would be fruitful to test in detail, as a large- 
scale disk warp may be excited by a planet on an inclined 
orbit (e.g., Lubow & Ogilvie 2001). 

We stress that our axisymmetric model assumes that 
the dust and gas are well mixed, it does not take into 
account the settling of the larger grains in the disk mid¬ 
plane. Moreover, having no observational constraints at 
longer wavelengths, we use conservatively a gas plus dust 
mixture that reproduces the images and SED of HH30 IRS. 
Other grain properties combined with dust settling might 
possibly account for some of the observed color gradient. 

In summary, while we can reproduce many of the fea¬ 
tures of the FS using a symmetric disk plus envelope, our 
models should be extended to include non-symmetric ge¬ 
ometries, specific dust properties, and dust settling, to ex¬ 
plain the shortcomings of the model we have presented. 
The PIGMap method introduced in this paper is a power 
tool to compare the predictions of these improved models 
with the strong spatial constraints brought by our obser¬ 
vations. A more extensive coverage of the SED of the ES 
from NIR to MIR and millimetric, combining ground and 
space observations (e.g., HST, VLT/NAGO-LGS, SIRTF, 
Plateau de Bure Interferometer), will further help to con¬ 
strain the dust properties (e.g., D’Alessio, Galvet, & Hart¬ 
mann 2001; Wood et al. 2002) and the disk structure of 
this object. 
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